The low energy spectrum of finite size metallic SWNTs 
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The electronic spectrum of metallic finite-size single-wall carbon nanotubes at low energies is 
derived. It is based on a tight-binding description for the interacting pz electrons. Not only the 
forward scattering parts of the Coulomb interaction, which are diagonalized by bosonization, are 
considered, but also all other processes becoming relevant for small diameter tubes. As a consequence 
of the substructure of the underlying lattice, a spin 1 triplet is found as ground state if the exchange 
splitting is larger than the branch mismatch, a spin singlet otherwise. Moreover the excitation 
spectrum is calculated. 

PACS numbers: 73.63.Fs, 71.10.Pm, 71.70.Gm, 73.23.Hk 



Single-walled carbon nanotubes (SWNTs) are one of 
the most prominent examples for the realization of ID 
electronic systems with orbital degeneracy in nature. A 
proper description of the low energy regime has to take 
into account the Coulomb interaction between the elec- 
trons and the corresponding correlation effects. Lut- 
tinger Liquid behaviour, leading to power-law depen- 
dence of various transport quantities, has been predicted 
theoretically for metallic SWNTs of infinite length 0| 
and was confirmed experimentally 0, 0|- Considering 
tubes of finite length, Kane et al. [5|| have derived the 
discrete energy spectrum of the collective spin and charge 
excitations and its dependence on the forward scattering 
part of the Coulomb interaction. Moreover, for SWNTs a 
charging energy of the order of the level spacing is found, 
leading to the observation of Coulomb blockade in SWNT 
quantum dot devices @, 0, S, S, Edl • The accompanying 
two- or fourfold periodicity of the Coulomb diamond size 
can be understood by including the spin and band degrees 
of freedom. A theory of transport through SWNT quan- 
tum dots incorporating the mentioned forward scattering 
interaction processes was worked out in [ll|. As discussed 
below, the restriction to forward scattering terms is justi- 
fied for large diameter tubes only, whereas the remaining 
interaction processes become more and more important 
for decreasing tube diameters. They lead to exchange 
effects and a modification of the excitation spectrum. 
Within a meanfield approach Oreg et al. predicted 
exchange effects favouring spin alignment. The size of 
the exchange splitting on the SWNT ground states was 
measured by recent experiments [3, Sill- 

A fundamental question regards the nature of the 
ground state and excited spectrum of correlated ID sys- 
tems. In a milestone theorem Lieb and Mattis |l3l | 
demonstrated that for a ID Hubbard model with near- 
est neighbour hopping, the ground state must have spin 
or 1/2. They left open the question for systems with 
orbital degeneracy. This is the case for SWNTs due to 
the substructure of the underlying honeycomb lattice. 

In this work we go beyond the mean field approach 
and derive the electronic structure of metallic SWNTs in 



the low energy regime from a microscopic model. Beside 
the long ranged forward scattering terms of the Coulomb 
interaction, which are exactly diagonalized by bosoniza- 
tion, we take into account all other processes becom- 
ing relevant for small diameter SWNTs away from half- 
filling. Because of scattering processes involving the or- 
bital degree of freedom, we predict a spin 1 triplet as 
ground state, if the exchange energy exceeds the en- 
ergy mismatch between the two electron branches and 
if 4m -I- 2 electrons occupy the SWNT. Moreover the ex- 
citation spectra are calculated. The large degeneracies of 
the discrete energy levels, as obtained by only including 
forward scattering terms, are partly lifted and the spectra 
become quasi-continuous when going to higher energies. 

The low energy Hamiltonian of metallic finite size 
SWNTs Without loss of generality we focus on arm- 
chair SWNTs As basis of our future discussion let 
us recall the minimal model Hamiltonian to describe 
finite size armchair SWNTs at low energies derived in 
[lH. It is based on electrons localized on the graphene 
honeycomb lattice which contains two carbon atoms, 
p = ±, per unit cell. Let the SWNT axis be along 
the X direction and ignore for the moment interactions. 
Then, by imposing periodic boundary conditions along 
the circumference and open ones along the tube length, 
eigenfunctions are standing waves (prK{r) with the 
branch or pseudo-spin index r = ±. The wave number 
K describes the slowly varying oscillations of iprKi'r)- 
The finite tube length L leads to the quantization 
condition k = 7r(mK + A)/L, € |A| < 1/2. 
The parameter A is responsible for a possible energy 
mismatch between the branches r = ±; its value depends 
on the length and on the type of the considered SWNT 
Explicitly, iprnir) can be decomposed into its 
contributions from the sublattices p = ±, ifrKi'r) — 
71 Ep=± fpr Ef=±k„ sgn(F)e-g-(^)«>pf (f), where 

f+r — l/v^ and f-r — The functions (pppir) 

describe fast oscillating Bloch waves at the two in- 
dependent Fermi points F = ±Ko of the honeycomb 
lattice, ippF{r) — N^^^'^ J^R^^^^^xi^ ^ R — ^p)i where 
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xir — R — fp) is the Pz orbital at lattice site R on 
sublattice p and Nl is the total number of lattice sites. 
The non-interacting Hamiltonian then reads 



Hq = hvF ^ 



where vp ~ 8.1-10^ m/s is the Fermi velocity of graphene 
and CrcTK annihilates an electron in the state \(prK,) 
Introducing the slowly varying ID operators iprarix) 
defined along the tube axis, 



1 ^e-gn(J^)K.^^^^^ (1) 



[/;/,—/,—/]. As we will show in the sequel the inter- 
action part V of the Hamiltonian is of the form V = 

Es.=/.f,,u T,SF=f,b Es„=/ Vs,s,.s. , where 

Vs^SfS^-^'^ X! J J dxdx'U[r][F]{x,x') 

We start with the Sr scattering types. From ^ it is 
evident, that the effective interaction potential is only 
nonzero for r2r3 = rir4. Hence we have to distinguish 
between the following two cases. 



and integrating over the coordinates perpendicular to the 
tube axis, the interaction part of the Hamiltonian be- 
comes effectively one dimensional. We find 



a) ri = r4, r2 = and b) ri 



V = 



\Y1 Y1 ti^S^i^i) I 1 dxdx'U[r][F][x,x') 

aa' {H,[F]}«=1 
X ^\aF^ Wl^l^a'F^ {x')^r3^'F3 {x')i^r^aF, (x) . (2) 



Here X]{[r] [f]} denotes the sum over all quadruples [r] — 
iri,r2,r3,r4) and [F] = {Fi, F2, F3, F4). Assuming that 
the wave functions (/5pj?(r) and Lp-p,p{r) do not overlap, 
the effective ID Coulomb interaction potential 



U[r][F\{x,x') 



1 r 
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C/[^'f™(x,x')(l + rir2r3r4) + 

C/[^f(a;,a:')(r2r3 + rir4)j , (3) 



can be separated into an interaction for electrons on the 
same (intra) and on different sublattices (inter), where 



C/™p/"'*''^(a:, x') j j dWd^r'j^ 

X ^*pFA^)'P*±pF2i^')'P±pF3{r")^pFAr)U{r- f'), 

and U{r — r') is the Coulomb potential. For the 
actual calculations we model U(r — r') by the so 
called Ohno potential. Measuring distances in units 
of A and energy in eV, it is given by U{r — r') = 



J7o/yi-t-(C/oe|?^-r'|/14.397) eV A reasonable 

choice is L/q = 15 eV Jl^. The dielectric constant is 
given by e w 1.4 - 2.4 [ij. 

The relevant scattering processes To proceed, it is 
convenient to introduce the notion of forward (/)-, back 
(6)-, and Umklapp (u)- scattering for an arbitrary in- 
dex quadruple [/] associated to the four electron op- 
erators in ([2]). With Si being the type of scatter- 
ing process for the quantity /, we denote a quadru- 
ple [/, ±7, ±/, /] by [I]sj^f±, while [/](, is equivalent to 
[/,—/,/,—/]. Finally Umklapp scattering means [/]„ = 



In case a) we have Sr — f and U[r]f[F]{x, x') is the sum 
of intra- and inter- lattice interaction. Case b) comprises 
Sr = b and Sr = w, with U[r]^^^[F]{x,x') being the differ- 
ence of the two types of sublattice interactions. Now we 
look at the allowed Sp processes. Although not dealing 
with an infinite system, after the integrations along the 
tube axis in ([2]), only terms with Fi + F2 — F3 — F4 — 0, 
i.e. the Sp = f and Sp = b processes, dominate as 
a consequence of the approximate conservation of quasi 
momentum. All other processes, including Sp — u, have 
very small amplitudes. Finally, regarding the spin in- 
dex, only Sa — f terms are allowed, since the Coulomb 
interaction is spin independent. 

Additionally, away from half filling, only terms with 
riFi + r2F2 — r^F^ — r^F^ = are relevant in ([2]), due 
to conservation of the quasi momentum, arising from the 
slow oscillations of the ID electron operators. 

Since C/[^j*™ and U'^^"'' differ only at the length scale 



of the lattice spacing the interaction potential for case 
b) [-^ Sr = b, u) is generally short ranged compared to 
the slowly varying electron operators tpra-F, whereas in 
case a) {-^ Sr = /), this is only true for the Sp — b 
terms. Hence for Sr = b, u or Sp = b the corresponding 
interactions become effectively local: 



Vs^SfS„ ■= Lu 



E 

{Hs,,[F]spJ<t]fs} 



X / dxi^laFii^)^l2a'FA^)^r3a'Fd^)A4'yFAx), (5) 



where the coupling constants are given by u^'' '^^ = 
1/ {2L^) J J dxdx'U[r]sr..[F]spix,x'). It holds m''^^ = 
y^uSp — . h'^Sf g^jjjj .^^g (define u+ :— u^^. The ratio 
y^Sr Sf j jg independent of L but scales like l/d where 
d is the tube diameter, such that these processes become 
neghgible for large diameter tubes. Numerically we find 
0.22 [0.28]f A and u^^ = 0.14 [0.22]f A for 



,Ah 



1.4 [2.4]. 
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Density-density interactions In the next step we in- a small perturbation to the Hamiltonian Hq + Vpp. We 
troduce the quantities Vpp and V^pp that collect all in- calculate the low energy spectrum and the correspond- 



teraction processes that are of density-density and non- 
density-density form, respectively, such that in total 
H — Ho -\- Vpp -\- Vnpp . Since the short ranged inter- 
actions are treated as local, we obtain 

^pp ^fff + ^/+&/+ + V'bf+/bf+, 

where the dominating part of the interaction away from 
half filling is the long -ranged term Vf f f. Using bosoniza- 
tion techniques [l7|, the Hamiltonian Hq -j- 
diagonalized and we find 

Ho + Vpp^J2J2 ^3^'}"']sq"'JSq + J^cA/'c 



Vpp can be 



jS q>0 



Mra (eo - '""'") + - T^^-ro 



(6) 



where we have defined J := 2{u^^ + u^''). The first 
term describes discrete excitations created/annihilated 
by the bosonic operators ajSq/alsq- ^o^'" chan- 

nels jS — c+,c—, s+, s-\- are related to total and rela- 
tive (with respect to the r index) charge and spin ex- 
citations. With the level spacing of the free electrons, 
Eq :— hvpj^, and e^q := e^nq, q = Uqj^ the relations 

ec+q = eogv^l + 8Wq/eo, e^/c-q = eo<;(l - u^''/eo) and 
Es+q — £0(?(1 + u^'^/^o) hold. Due to the dominating 
Vf f f contribution, the ratio gq :— eoq/ec+q is strongly 
reduced for small q {gq « 0.2) and for large q, gq tends to 
1 Small corrections due to the coupling constants 

and u"^^ have been omitted. Finally, 



Wq = -^ dx 

^ Jo Jo 



dx' cos{qx) cos{qx' )U[r] fiF] f 1 x')- 



The remaining terms in ([6]) are fermionic contributions, 
accounting for the energy cost of changing the number 
of electrons in the different branches (rcr). The opera- 
tors Nra count the electrons in (rcr) and Nc. — J^ra^ra- 
The single summands account for (in the order of ap- 
pearance) the Coulomb charging energy Ec — Wq, the 
shell filling energy (because of Pauli's principle), a possi- 
ble energy mismatch £a = sgn(A)eo min(2 | A| , 2 | A| — 
1) between the r branches if |A| ^ 0,1/2, and a 
favourable spin alignment of electrons with different r 
due to Vft f+ /b f+ . Note also that the shell filling energy 
is modified by the attractive contribution — u+ due to 
Vf+i,f+. The eigenstates of Hq + Vpp are \N,rn) :— 

UjSq[a]sq) /^/^■\N,0), where \N,0) has no 

bosonic excitation and N = (Af_|, 7V_|, iV_|_|, A^+j^) de- 
fines the number of electrons in each of the branches (ra) . 

Non- density- density processes In the following we 
concentrate on the situation away from half-filling, where 
the non-density-density processes Vnpp can be treated as 



ing eigenstates of the full Hamiltonian Hq 4- Vpp -I- Vnpp 
by evaluating the matrix elements ( N,rn\Vnpp\N\rn' ) 
and by truncating the Hilbert space at a certain eigenen- 
ergy of Hq -\- Vpp. Near half-filling the strength of Vnpp is 
highly enhanced and the truncation procedure therefore 
questionable. We obtain 



Vnpp - Vf+bf- +Vbf+/bf- +Vuf-/bf- 



(7) 



Low energy spectrum Our truncation scheme to find 
the low energy spectrum is to only retain the energetically 
lowest lying states of Hq -f Vpp which have no bosonic 
excitations. Using ^ and ([5]), we get for the matrix 
elements of the contributions on the r.h.s of ([7]), 



[N,0\Vs^s,sjN',Oj^-u'''-S- 



E 



X 5- 



J:^^! sgn(Fi)Ki~Y.U3 sgn{Fi}Ki,Oi 



(8) 



where the Kronecker-(5 results from the integration along 
the tube axis in ([5]). Note that the states |iV,0) are 
eigenstates of Hq -\- Vpp and not of Hq alone. Hence the 
evaluation of |[8| in general is not straightforward. But 
for the processes relevant away from half-filling we get 
the same result and physical insight if we consider for 
the calculation of ^ \N, ) as the Fermi sea filled up 
with Nra- noninteracting electrons in the branch ra, i. e. 
as an eigenstate of Hq , as we will discuss elsewhere . 

In the following we focus on the case Nc — 4m -I- 2. 
As truncated basis we use the states \N,0) with N = 
{m -\- l,m -\- l,m,m) and permutations. In the following 
we denote \{m -\- l,m + l,m, m), 0) by |ti, — ) etc. Using 
([6]) and ([8]) the interacting Hamiltonian, restricted to the 
states |T,T), IT,i), li,T), in,-) and |-,n), is 

H = EQ^4m+2 + 



I 





r 2 


2 " 



V 

with i?0.4m+2 



- eA 
J 

2 



(9) 



£a / 



+ (2m2 + 2to -I- 1) (eq - u+) - 
J(m} -t- m) -I- 2u+TO. Diagonalizing H, we 

find that its eigenstates are given by the spin 
1 triplet |T,T>, U,i>, l/V2(|T,i> + U,T)), the 
spin singlet l/%/2 (|t, i> - U, T)) and the two 
states (c2/2 -f l)-i/2(ci/2|n,-)±|-,n)), where 

Ci/2 = 4/ 



( J/2)2 ± £a ) • Relatively to i?o,4m+2, 
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the eigenenergies are —J/2 for the triplet states, J/2 
for the singlet state and ± -y/e^ + (J/2)^ for the 
remaining two states. Thus, under the condition 
J/2 > + (J/2)2 - M+, i.e. for a small band 

mismatch ea ^ J/2, the ground state is threefold de- 
generate and formed by the spin 1 triplet, otherwise by 
{cl + l)-^/"^ {c2 ITi,-) - |-,Ti))- The low energy spectra 
for the two cases ea — and £a ^ J/2 are shown in 
Fig. [T] for a (6,6) armchair nanotube (corresponding 
to c? = 0.8 nm) with an assumed dielectric constant 
of e = 1.4. The obtained values of J = 0.09eo and 
= 0.03^0 are in good agreement with the experiments 
where nanotubes with £a ^ ^ were considered. 
In accordance with the discussion above, Moriyama et 
al. [l3| could identify the ground state to Nc = 4m + 2 
as a spin singlet and also a spin 1 triplet was found. 
Not observed so far has been the spin 1 triplet as ground 
state and the mixing of |ti, — ) and | — , ti), as predicted 
by our calculation for nanotubes with ea ^ J/2. We 
emphasize that all the exchange splittings here result 
from non-forward scattering processes with respect to 
the band index r. In the considerations of Lieb and 
Mattis [3| such an additional pseudo-spin degree of 
freedom is missing and so we conjecture that this is 
the reason why their theorem can not be applied in 
our case. The meanfield result of Oreg et al. [l3| for 
the ground state structure to Nc — + 2 essentially 
can be obtained by setting the off-diagonal elements 
in |[9l to zero. Therefore if £a ^ -j, the meanfield 
approach yields the same ground state spectrum as 
our work, but with different degeneracies: Instead of 
having a threefold degeneracy of the spin 1 triplet 
and no degeneracy for the spin singlet, a twofold 
degeneracy of the states |t,T) , and |T,i), |i,T) 

respectively, is obtained. Additionally for ea < J/2, the 
meanfield theory is not capable of predicting the mixing 
of the states ITi,— ) and | — ,Ti) with the accompanying 
exchange splitting. The ground states of Hq + Vpp to the 
charge states Nc ^ 4r7i + 2 do not mix via Kipp and thus 
the corresponding energies can be determined easily by 
using eqs. {6]) and ^ ior N = N' . 

Excitation spectrum So far no bosonic excitations 
were involved in our ground state examination. In order 
to discuss the excitation spectrum oiH — Ho+Vpp-\-Vnpp, 
the matrix elements of ( N ,rh\Vnpp\N' ,Tn' ) must be de- 
termined. The detailed calculation, based on bosoniza- 
ton techniques, will be given in a longer article After 
truncating the Hilbert space at sufficiently high energies 
and diagonalizing H we obtain for Nc — 4m -I- 2 the spec- 
trum as shown in Fig. [2l For comparison we also show 
the spectrum of the "standard" theory 0, [ll| as ob- 
tained by only retaining the dominating forward scatter- 
ing processes of Vpp. Most striking is the partial lifting 
of the huge degeneracies and the formation of a quasi 
continuum at higher energies. 

In conclusion we have derived a microscopic low en- 
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Figure 1: Low energy spectrum of a (6,6) SWNT for Nc = 
4m + 2. (a) For £a = the ground state is formed by the 
spin 1 triplet (^ ©). The states |Ti,— ) and | — ,Ti) mix. (b) 
For EA 3> J/2 the ground state is given by the spin state 
jtl, — ) . The spin singlet is indicated by ®. The interaction 
parameters are J — 0.09eo, u"*" — O.OSeo. 
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Figure 2: Excitation spectrum of a (6.6) SWNT including 
only forward scattering processes (— > x) and for the total 
Hamiltonian H = Ho + Vpp + Vnpp (^ <)). Note the logarith- 
mic scale of the a;-axis. The energy of the lowest c-|- excitation 
is 4.3 eo- All other interaction parameters are as for Fig. [T] 

ergy theory for finite size metallic SWNTs away from 
half fining, including non-density-density interaction pro- 
cesses, which become relevant for small diameter tubes. 
The ground state and excitation spectra have been deter- 
mined. In particular, we predict a spin 1 triplet as ground 
state for Nc — AiJi+2, if the energy mismatch ea between 
the different pseudo-spin branches is much smaller than 
the exchange energy J, a spin singlet otherwise. For 
£a S J/2 we furthermore find that the pseudo-spin is 
not conserved and the corresponding degeneracy is lifted. 
We notice that in an energy mismatch ea^ ^ was 
used to fit the data and, in agreement with our theory, 
a singlet ground state has been inferred from magnetic 
field measurements. Observation of the triplet ground 
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state is within experimental reach. 
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